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We study a particle on a ring in presence of various dissipative environments. We develop and 
solve a variational scheme assuming low frequency dominance. Our solution produces a renormal- 
ization group (RG) transformation to all orders in the inverse dissipation strength, and in particular 
reproduces known two loop results. Our RG leads to a weak dissipation parameter, for which a weak 
coupling expansion for the position correlation function shows a decay in imaginary time. 

o 

_ PACS numbers: 

^ ' I. INTRODUCTION 

fSJ . The problem of interference and dephasing in presence of dissipative environments is of significance for a variety 
of experimental systems and a fundamental theoretical issue. The experimental systems include mesoscopic rings 
embedded on various surfaces where Aharonov-Bohm (AB) oscillations can be measured^'^, and the related problem 
' ^ ■ of decoherence at low temperatures^. A different type of experimental systems are cold atom traps created by atom 
' chips^^^. The atom chip that produces a magnetic or electric trap for the cold atoms necessarily also produces noise. 
^ Our problem is then relevant for evaluating the interference amplitude of the cold atoms in presence of such noise. 
Qj As an efficient tool for monitoring the effect of the environment we follow a suggestion by Guinea^ to find the 

AB oscillation amplitude as function of the radius R of the ring; for free particles of mass M this amplitude is the 
. mean level spacing ~ l/MB?. Two types of environments were suggested to lead to an anomalous suppression, i.e. 
"j^ , a stronger decrease of the oscillation amplitude than a Caldeira Legget (CL) bath as well as a charge - metal 

' (CM) system, i.e. a charge on the ring interacting with a dirty metal environment. The CL system is of further 
I \ interest since it can be mapped to the Coulomb blockade problem^ii^ as well as to quantum dots at a distance from 
' metallic gate&iS. The Coulomb box problem is of further recent interest in view of data on the quantization of the 
Ch \ charge relaxation resistance^iii^ and related theoretical developments^^"— . 

O . The CL system has been extensively investigated by instanton method a^^d^ , by RG methods^'®, by a boundary 
I field theoryi^ and by Monte Carlo (MC) methods^i^ii^. All methods show that the effective mass, defined as of 

' . . ..... . 2 . . 

the particle increases exponentially with the dissipation strength a, i.e. B a^e^ with differences in the exponent 
jjL. In 2nd order renormalization group (RG) = —1— while instanton methods give either— jj. = —2 or— /x = —3; the 
boundary field theory with MC gives /i = —2. A variational approach-- indicated a nonperturbative regime at strong 
a. Since a — ^B? where 7 is a friction coefficient, a length scale tt/ ^/^ is identified''; this scale is a candidate for a 
dephasing length. 

The CM system was investigated by RG methods^ finding B ^ with /x' < 1 nonuniversal, while MC data^i 

shows /i' ss 1.8. Further MC simulations show that in fact ji' — 0, at least for weak coupling^. We study also a 
dipole-metal (DM) system, i.e. an electric dipole on a ring coupled to a dirty metal environment. This system can 
\ be realized by experiments on cold Rydberg atoms^. 
^-H . In the present work, extending our previous reporlj^l we solve these systems by a variational method, assuming 
low frequency dominance. We find that the variational method defines an RG scheme to all orders, reproducing a 
known RG equation^ to two loops in the CL system. In the CM and DM systems, for either a charge or a dipole, we 
find that the effective mass remains B/B? ~ i?" for large R, as for free particles. Our RG leads to a weak coupling 
dissipation parameter. The resulting action yields a weak coupling expansion for the position correlation function, 
showing a decay in imaginary time. This decay is generic to all finite B systems and indicates dephasing of an 
excited state. In the limit i? — > 00 the correlation probes degenarate states, however, the position correlation function 
does not decay in this limit, i.e. no dephasing. 

In section II we present the models. In section III we define our variational method and show that the effective mass 
B of the m = sector determines the curvature d"^ Eq/ d4>l.\o where Eq is the ground state energy and (j)x is the flux 
through the ring; this curvature is a measure of the Aharonov Bohm oscillation amplitude. In section IV we simplify 
the variational equation by assuming low frequency dominance, or equivalently logarithmic dominance. In section V 
we show that this method is equivalent to an RG scheme, and in particular reproduces the known RG equation to 
2nd order in the CL system. In general the variational equation contains terms to all orders and is therefore expected 
to be superior to a 2nd order RG expansion. In section VI we present explicit solutions for the CL and DM systems, 
as well as for a general case. Finally in section VII we study the weak coupling expansion showing a decay for 
the position correlation function. 
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II. THE MODEL 



In this section we derive the effective action in presence of a dissipative environment in terms of the the angle 9m (t) 
where r is an imaginary time. The index m specifies the winding number so that 

9m{T)=e{T)+2TnnT/l3 (1) 

where 0{Q) = 9{j3) has periodic boundary condition and /3 is the inverse temperature (/3 — > oo below). In presence of 
an external flux (px the partition sum has the form 

Z = ^ /"255/g27r»m0,-Si{e„.}-s„t{e,„} ^2) 

As shown by Guinea^, the form of such an action in presence of a general dissipative bath the effective action can be 
written in terms of a kernel K[0{t) — 0{t')] that is periodic and allows in general a Fourier expansion 

Si{0} = / dT^i — f 







rP rP 
Jo Jo 



sin^[7rr(T - r')] 



= o v r r -.(.01/2) 

and a„ depends on the type of bath. At r — 7> r' (or at high frequencies w) one can expand the sin^(...) in ([3]) and 
then Sint — >■ Oi^n'^n'n'^ j duj\uj\\9miy)\'^ , identifying a dissipative system. 

We consider now 3 types of environments and identify the coefficients a„. First is the Caldeira Legget (CL) 
environment. It has harmonic oscillators coupled linearly to the particle's coordinate. The effective action is well 
knowEk^^ for nonconfined coordinate R(t) 

™* 'i i sin2[^T(r-T')] 

where 7 is the dissipation parameter. When the particle is confined to a ring R(r) = R[cos 6'(t), sin 0(t)] the action 
becomes of the form of Eq. ([3]) with a single coefficient oi = 1 and a — ^E? . 

Consider next the charge-metal (CM) environment. It consists of a dirty metal that is characterized by its conduc- 
tivity a and diffusion constant D. The particle on the ring has a charge e and responds to the Coulomb potential of 
the metal V(R,{t)^ t). The metal is assumed to be a Gaussian environment, so that the interaction term (in imaginary 
time) of the partition sum can be averaged to obtain^! 



and with / = / (fq/{2TTf 



S„,t = ie^ (F(R(r),r)V(R(r'),r')) = ^e^T^ / ■ 

n •'1 



^ g^q■(R(r)-R(r'))-ic^„(r-r') j^g^ 

_q^e{i\u:n\,q) 

where the propagator of the scalar potential^^ is given in terms of the dielectric function e(j|a;„|, q) with w„ — 2iTnT 
the Matsubara frequencies. At low frequencies and momenta e(a;, g) = 1 + , valid at g < 1/^ where t is the 

mean free path. Hence l/e(j|w„|, q) sa (l^^iil + Dq^) / Ana] the Dq^ term yields an R(t) independent constant while 

sm [tt i t\ 

hence with kp the Fermi wavevector and r = R/ i the charge coupled to a dirty metal has 
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if(z) = 1- [4r2sin2(z/2) + l]-i/2 (g) 
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For r ^ 1, Qn ^ ln(r/n) for 1 < n < r and a„ for n > r. This model reduces to the CL one at r <C 1, where 
Ui single a„ survives. 

A 3rd realization of the action corresponds to the dipole-metal (DM) environment. Consider a particle with 
an electric dipole, whose direction is perpendicular to the ring, interacting with a metal. For the electric field 
Ez = dzV — ^dtAz, the propagator involves^^ [e{i\uJn\,q)uj^ + Dq'^]~^, which for q ^ can be expanded in uj^, 
hence it has no dissipative term ~ |wn|; we keep then just the dzV term. The interaction with the fiuctuating electric 
field E(r, r) is p i?^(R(r), r)(ir. A Gaussian average on the metallic environment then yields 

S^nt = hp' I I (a.y(R(T),r)a.V(R(r'),r')) = Ip'T^ / -^-i^^^e^^' (9) 



q q'^<^{i\^n\,q) 



Therefore 



K{z) = l-(4r2sin2| + l)-^/^ (10) 

Hence, for large r, a„ ~ ^(1 — ^) for n < r and a„ sa otherwise. Finally we note that a topological flux (p^ can be 
realized for an electric dipole^i. 

III. VARIATIONAL METHOD 

The partition sum can be rewritten by using the Poisson sum J2m di''^) ^ J^k I di'^) exp(27riif 0)(i0 so that 

/"-^ r 2 2 2 

d(l) / X'6>e^"'^<^+'^'°^~ '" T -Sl{9(r)}-S.„t{e(r)+2^0r//3} 

^ r dcpe^^'^^''+^-^-'^^^ (11) 
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The variational method for finds the best Gaussian approximation, i.e. 

^" = ^Eg"'k)i^k)P (12) 

SO that the free energy in = e~^"^ has the variational form 

fiF^ar = m + {S- So)o = ^ lnG(c^„) + [MR^Lul - G-\ujn)]G{ujr,)} + {S,nt)o (13) 

where (...)o is an average with respect to exp(— S'o) and Fq is the free energy corresponding to So- Since 

2(sin2{n[0(T) - e{T') + 27r0(T - t')//3]/2})o 

= 1 - cos[27rn0(T - t')/I3] exp{-n2([(6)(T) - 0{t')^])o/2} 

= 1 - cos[2nn<j){T - t')//3] expj-n^T^ G(w„)[1 - cos(w„t)]} (14) 



the interaction term becomes 

(5,„t)o =/3a Va„ ^^{1 - cos(2^n0r//3)e-"'^^" ^^'^^Ii-^^^"^)!} . 

The variational equation dFyar/SG{ujn) = is then 

G-\u:) = MRW + 2a V a,,n^ f ^^ I-cos(c^t) ^os(27rn0T//3)e-"' /(d"i/2-)G(^,)[i-cos(^ir)] 

"'0 



(15) 



(16) 
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When the hmit /3 — >■ cx) is taken a cutoff Wc may be introduced to control the short time behavior so that the r integral 
becomes . This cutoff represents a high frequency limit of the bath degrees of freedom. Alternatively, the mass 

term serves also as a cutoff since it leads to convergence of the duti integral in the exponent of 

In the following we will study the variational equation with (j) — 0. To justify this, we show now that the effective 
mass B of the cj) = system is indeed what is needed to find the Aharonov-Bohm oscillation amplitude at T — >■ 0. 
The effective mass is defined by G~^{uj) ~ Blu^ in the limit oj and is identified from Eq. ([T5]) at /3 — > oo as 

1 roo 

B = MR^ + -aY,^nn^ ^^g-nV(<i-/2x)GM[l-co.(-r)] (^7) 

n 

The form (fTTj) implies that the fiuctuations < 0^ >^ /3, hence the factor cos(27rn0T//3) — > l + 0(l//3) in Eq. (fT6|) and 
the effective mass B is 4> independent. It is also necessary to check that the r integrals converge: indeed at r — > oo 

/ dujG{uj){l - cos(wt) w tM — + / dujG(uj) - t/B 
Jo Jo 2i? Ji/^ 

hence a factor e^" '^^^ assures the convergence of the t integrals. 

The Aharonov-Bohm oscillation amplitude is usually measured^i^ by the curvature of the free energy at (j>x — 0; 
since at (/> = we have dG{uj)/d(f> — (from parity in <j), see (fT6|) . and from analyticity in <j>) and dF^ar/dG = (the 
variational condition) we obtain from Eqs. (|13|15|17p 



''^^-lo^m^\o = '-^iB-MR^). (18) 



, , , - I n — - i\, I ^ ^ 



The effect of in the partition sum Eq. (fTTj) is therefore to replace the factor M (jp' / j3 by 27r^_B0^//3, i.e. the 
response to an external flux is that of a free particle with a mass renormalized to B. Higher order terms produce only 
subdominant behavior in 1//3, e.g. one expects a 0*//3'^ term. Our task is therefore to study the = system and 
find this renormalized mass. 

IV. VARIATIONAL EQUATION 

Before studying the full equation, it is instructive to study its perturbative regime. The lowest order is obtained 
by neglecting the exponent in Eq. (|16|) . leading to 

G"^(a;) = A/i?V +7ra;a^a„n2 uj<uJc (19) 

n 

This identifies the cutoff uJc below which dissipative term dominates, 

- MR^ 

Consider next lo <^ lOc but still ln(ti;c/w) < 1. The next order in perturbation is obtained by using Eq. (|19p in the 
exponent in Eq. (jl6p and expanding this exponent, 

G~^(w) = TTCjay^ a„n^[l — =^ :r In — 1 = Trwa a„n^[l In — 1 ln(a;c/a;) < 1 (21) 

a > cimTO UJ ^ an lj 

/ ^ III (VI 



(22) 



where the mass term is ignored for w < and k is a geometric parameter defined by 

The sums in p2|) can be evaluated for each model from the 2nd and 4th derivatives of if (z) at z = 0, leading to 

= ^ CM 

= (23) 
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A significant perturbative regime is possible for aK ^ 1. This strong dissipation condition can apply to the CL 
model if R is large, though one needs to make sure that the CL model is still valid in that case. For the CM or DM 
models k is bounded by a number ~ 1 so that a S> 1 is needed. For usual dirty metals kp£ 1 so that for charge 
coupling with Eq. the condition is not satisfied, unless the particle on the ring has a charge e* ^ e. On the other 
hand, the dipole case may have a large a in Eq. ([T0|) for large dipoles, e.g. in Rydberg atoms. In the following we 
use the form (|21|) as a boundary condition for the full variational solution. 

We proceed now to variational equation, that includes the significant range of <C Wc- It is convenient to study a 
derivative of Eq. (ITCl) 

du} n "'0 

The bare mass M serves to define ujc and then the Mi?^w^ term in Eq. (|16p is neglected at w < ojc- If uj is sufficiently 
small then sin(a;r) can be expanded leading to a ^ cj term. Wc therefore assume the form 

G^^ = /(w) (jjQ < uj < ujc 

Q-'^ = Buj^ uj <u}o (25) 

The solution for f{uj) needs to satisfy boundary conditions, whose a dependence is determined by the perturabtive 
expansion Eq. (HI]), 



/'(We) =™E„a«n'-??(a) Ti{a) = 1 + — + 0{1/ a") 

n^c) =7rE„a„n2.^ C(a) = i + 0(l/a) (26) 

We proceed to simplify Eq. (|24p . For w > wq the oscillating sin(wT) in Eq. leads to a cutoff t < l/rjiLo, to be 
determined by matching to the perturbative regime. Hence 

/'(w) = 2wa Va„n2 r^"" drg-"' /o^^('^"i/2.)«("i)[i-™^("i^)l (27) 

The range involves G~^{lu) — AIR^ui^ and contributes ^ 1/{MR'^uJc) = [^ct X^n '^""■^l^^ which is neglected for 
a ^ 1. The r integration is dominated by t l/r/iw, hence 1 — cosojit 1 — cos(a;i/77iaj) is replaced by ijj\/2uP'rj\ 
for wi < a; and by 1 for uji > lj. This rough separation is to be justified by our main assumption that " duji/ f{uji) 
dominates this integral due to the low frequency decrease of /(i^i). The terms from ujq < uji < uj, as well as those 
from LUi < wqi can be neglected if 

/ condition (i) . (28) 



Note that the 2nd term on the left near cjq is ~ I/Bcjq, while near ojc it is ^ l/a and negligible for large a. We 
are interested in nonperturbative contributions, i.e. the range Inwc/"^ ^ 1 and in particular at w = wq- In terms of 
a;2 = I/t we obtain 

/ (w) = 2wa > a„n / — ^-e ■'"a'^s ^' (29) 



where as above, the precise location of the ri2Uj cutoff should not be significant. The lu2 integration is dominated 
by its lower cutoff r]iuj so we expect that the exponent can be taken out of the integration with the replacement 
i]2UJ2 — >■ r]ir]2U}. More precisely, taking a derivative of ([^ leads to 

/'(w) = TTfi{a)a ^ a„n2g-n^ J-" du^/^f(u^) _^ ^/"(^) (30) 

n 

and rji — 2/(7r7y) and 7717^2 = 1 are chosen. The coefficient f]{a) is to be determined by the boundary conditions (|26p . 
To further simply the equation we assume now 

/"(w) < condition (n) . (31) 
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leading to our main equation for f{uj), 

/'(w) = 7r77(a)a ^ a„n2e-"' (32) 

n 

The coefficient here is i]{a), consistent with ([26l) . Below we actually find that condition (ii) is not always satisfied, 
and then we return to solve Eq. (pO)) instead of (|32|) . 



Finally, consider w < cjq- Eq. ([M]) has then on the left -£^G ^{oj) = 2Bu! while on the right it has the requested 
^ cj form, except for a term where 

n V- 2 /•'/''^" , , 1- costal r , dui 
Yi = 2oja } a„n / o(Texp[ (/ duji + / — — -)J (33) 

Since r > I/ojq dominates, — coswir)c?a;i hence 

h = 4c^Ba^a„(e-"'/2B^° - e-»V2B.i-)g-"^ ''"i (34) 

n 

The essential singularity in uj is negligible for < when 

BujQ > 1 condition (iii) (35) 

The remaining term at w < cjq identifies B and leads to a matching condition of the form /'(wq) = tf Bu^. Continuity 
of derivatives yeilds 77' = 2, though we expect that the precise value of 77' will not be significant. 

This completes the derivation of the equations for B and /(w). Eq. ([50]) or ([5^ are to be solved with the boundary 
conditions in (in case of p2p only the first two conditions are needed). Furthermore, the matching conditions at 
LjQ are 

/(wo) = Bui 

f'iiuo) = v'Bluo 77' « 2 (36) 



V. RG PROCEDURE 



We present here an approximate solution of the variational equations by an RG method, which in some case (the 
CL case, see below) should be very close to exact. The idea is that an uj < lUc can serve as a new cutoff provided that 
the coupling a is renormalized into a{uj). The boundary conditions ()26p become therefore 

f{u) = TTUja{uj) a„n^ 

n 

f'iuj) = -Kajuj) a„n^ ■ ??[a(i:j)] 

n 

/ (w) = TT > a„n • (37) 

n 

The number of needed equations depends on the order of the differential equation for f(j-o), e.g. for Eq. (1321) only 
the first two equations in (j37p are needed. The functions 77(a), C(a) are known as an expansion in 1/a. As we find 
below, these functions can be determined explicitly by the variational equations. 
Taking a derivative of the 1st equation in p7l) yields a recursion relation for a{uj), 

^ aiuMaiio)) - 1] (38) 

Hence the boundary condition function i](a) determines the flow of the renormalized a{uj), i.e. it generates the RG 
flow to all orders for which 77(a) is known. In particular the flow terminates when 77(ac) — 1, i.e. ac is a fixed point. 

Before proceeding to solve for 77(a), we show that the RG is equivalent to a solution of the form /(w) = ujg{K{a)uj), 
so that all the a dependence is included in the function K{a), i.e. the function g(x) itself is a independent. This 
property is exact for our variational equation for the CL system (see below and Appendix |B]) . For other systems the 
scaling function needs to be identified separately, as e.g. done in section VIC for the CM system. 
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The first boundary condition from ([26)) is g{K{a)u!c) = t^ol 'Ylm '^nn^, fience K' {a)uJcg' {K {a)ujc) = tt o-nn^ ■ Tlie 
second boundary condition is then 

/'(wc) = g{K{a)ujc) + K{a)ucg{K{a)ujc) = {a + ^jf^)'^ """^ = 7r«??(«) 

^ ^(a) = l + -|7rT (39) 

In g{K{oi)uj) one can vary either a ov uj with identical effects if K{a)ijj — K{a(uj))LOc (see also Appendix B) which 
by d/duj yields 

K'[a{uj)]^^^^^u^ = K{a) = K[a{uj)]— (40) 
aw w 

and with p9p the flow (1551) is reproduced. We note also that the scaling functions ri{a), C{a) ([57)) can also be 
determined by rewriting the differential equation for g{x) in terms of y{g) — x{g)g' {x{g)) and its derivatives. This 
is possible under fairly general conditions, e.g. that 5(2;) is monotonic and that the differential equation for g{x) is 
homogenous (i.e. contains only x"g^^\x) terms). 

So far the RG flow was determined in general, without a necessity to identify the relevant differential equation. We 
proceed to show that the differential equation for /(w) determines the function 77(a) completely and therefore also 
the flow of a{uj). Consider Eq. ([5^ that leads to 

f (cjc = 7r?7Q " = TTrj—!^ 41) 

TT]-(UJc) KUJc 



Differentiation of the second equation in ([37| leads to 

where (arj)' — -^[a'r]{a)]. Equating the last equation at cUc with ([4T|l leads, in terms of £ — Incj, to 

'iU--^----^ + Oil/a^) (43) 
at K[ariy k n^a 




1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — I — |— 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

l/oK 



FIG. 1: Solution of Eq. (|44p for the scaling function rj as function of na. The branch at 77 < 1 is not accessible with the initial 
values of (|26|l . Note that k = tt^ for the CL system. 
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To obtain the expansion we use the perturbative form of rj in ((26)) . Remarkably, the result p3)) is precisely the two 



loop RG result For the CL system-^ (with g = in the notations or Ref^ and k 

that the same perturbative form in Eq. (|38| yields only the first term 1/k. 
Comparing Eqs. p8l43p yields 



TT for the CL system). Note 



r, = l 



Kat] 



1 



2 dri 



(44) 



This relation generates a large a expansion with the leading form ry = 1 + {Ka)~^ + 0{a)~^, consistent with the 
perturbation expansion Eq. (1261) . It is remarkable that the perturbation expansion allows for an asymptotic expansion 
of (l44l) . i.e. a different form of ?y(a) in ((26)) would not allow such an expansion. 

Fig. 1 shows the solution of this equation, with the exact analytic solution given in appendix A. Note the turning 
point at ?7 = 1, 1/(q:k) — 0.742. This corresponds to a fixed point at ac, i.e. if this point is reached at a frequency Wa 
then at io < uja a{uj) = ac remains constant and f{io) — UcU). This behavior is in fact inconsistent with the assumed 
form ([24)) . Another difficulty is that continuity of f {wq) / /{uq) — 77[a(a;o)] = rj' needs 77[a(wo)] ~ 2, which is not 
achieved in Fig. 1. 

In the next section we evaluate /(w) itself and show that the solution based on Eq. (|5^ does not satisfy criterion 
(ii) below some low frequency uj^ > uja- In the latter range one needs to address Eq. pop . we note that the Lof"{uj) 
term in (|30)) is small at the initial range of to, e.g. at ujc it is 0(l/a) relative to the /'(w) term. Therefore, to be 
consistent with the terms neglected due to the criteria (i), we need to start with Eq. and only at the frequency 

~ LOb we shift to Eq. ([50)). 

We proceed to study the RG form of ([^0)) . Taking a derivative of the second equation in ([57)) at cj = cjc yields 



C(a) = a[77(a) - l](ar;)' 



(45) 



Eq. (|30)) taken at ojc yields fj{a) = 77(a) — C{a)/a. Next we evaluate /'"(wc) in two ways: First, by taking a derivative 
of Eq. ([50)) that leads to /'"(wc) = — 7r77(Q;) a„n^/ (kw^). Second, by taking a derivative of the third equation in 
([57)) . Equating these two forms leads to, finally. 



77(a) 



C(a) 



= K{C(a) - C"(a)a[?7(a) - 1]} 



(46) 



Together with ()45|) this is a 2nd order differential equation for 77(a). We solve this equation by matching at some a 
to the solution of ([^H) . as shown in Fig. 2. Curiously, ()46p has an exact solution 7; = 1 + l/{Ka) which gives the one 
loop solution in ()43l) . As mentioned above, we apply this solution only below some low frequency ~ wq, to be studied 
in the next section i.e. 77 = 1 + — , does not have then the proper boundary conditions. 




FIG. 2: Solution of Eq. (|46p . upper line, for the scaling function 77 as function of Ka. The initial values for this solution are 
taken from a point on the solution of Eq. (1441) . lower line. 



VI. SOLUTIONS FOR VARIOUS SYSTEMS 



We present now explicit solutions for /(w) and study the validity criteria. We start with the mathematically 
simplest case, the CL system. 



9 



A. Caldeira Legget system 
Considering Eq. ([32)1 we obtain by differentiating 



which upon integration yields 

f{u:)=TT-'HKf{Lo)] (48) 
where K is an integration constant. Further integration yields 

li[Kfiij)] = Tl-^K{iU - iUa) + ll[Kf{^a)] (49) 

where at Wa, a{uj) reaches the fixed point of Eq. (j44p i.e. a(wa) = etc anticipating that this equation is not valid all 
the way to uja\ here li{x) — is the Log Integral function. K is determined by the ujc values 



so that 



and from /{lOc) = Traujc we have 



fiiUc) = ??7ra = TT-i HKfiiu,)] (50) 
Kfico,) = e"-'" » 1 (51) 



K = . (52) 

naujc 

Eq. (j48l) a.t Lj — oja yields 

Kf{u;a) = e''''^^"" (53) 
For B and iOa we need to solve the coupled equations 

dz K 

= — (Wc-Wa) 

tBuia inz TT 

KBloI = e^'''^-^"". (54) 

An explicit solution requires an asymptotic expansion of li{x), which is provided by our RG method. As discussed 
in section V, the solution has the form /(w) = Lug{Kuj) such that g{x) is a independent. Eqs. (|40I52[) then yield/(a;) = 
T:uja{uj) where a{uj) the solution of 

7r^a(w)jj[a(w)] 

Ku; = — . (55) 

7Ta(U!) 

Inverting this relation we find 

f{Lo) = T:ua{Lo) = — \n{na{uj)Kuj) = — ln(— \n{7ra{uj)KLo)) = — ln(— ln( — ...)) (56) 
TTT] nrj rjTT nrj rjir rjir 

and at least two In embeddings are needed for a large a solution, i.e. 

/(.) = i^ln[^ln^] + 0(-^). (57) 

The boundary condition at uja is Kua — e'' ^^'^•^ /Buja so that g{Kuja) — Buja becomes 

^rj TrBuJa 

ff(-g^)=SWa. (58) 
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This equation does not involve the large parameter a, hence Buja ~ 1, KuJa ~ 1, and the effective mass at scale LOa is 

B^-!-^- . (59) 

UJa aUJc 

We note that Eq. ([SS]) implies that a{u!a) = 0{1) i.e. in the vicinity of the fixed point ac = 0.14. 

We check now the conditions (i)-(iii) for lo near uJa- For condition (i) we use Eq. p2|) at a; = LUa, where f'{uja) = 
TTr][a{uia)]a{u}a) = 0(1), hence 

- In a (60) 



7r/(wi) 

Hence the condition (i) is satisfied only for In a ^ 1. The condition (ii) corresponds to a;/"(cj) — u f (uj) / tt f (u) ^ 
/'(oj), hence 7r/(aj)/a; = 7r^a(a;) ^ 1. At this condition fails. 

We consider therefore the previous solution as valid only down to a frequency to be determined below. At 
a; < Wb we use the more complete Eq. (30). As seen in Fig. 2, below LUb the slope ri{a) increases rapidly towards 
the value rj' 2 which determines ujq. A numerical fit to Fig. 2 and use of (38) yields a weak a dependence, i.e. 
Wo ~ tJb/o:^ with X ~ 0.2. Neglecting this effect, we identify ujq « uib and check the various conditions. Consider first 

duji f'{(^b) , a 
— -Trm = 7rin' 



/(wi) na a{uJb)r][a{uJb)] 



Since Bojq w f{u!b)/ujb — iraiiOb) condition (i) is satisfied for any choice of Wb such that a{ijJb) ^ 1, e.g. a{ijJb) = ol'^ 
with v < 1. 

a{LOb) 

with ri,r](a{ujb)) ~ 1; since < 1 a{ujb) < « and loq <C uJc and provides a huge range where Eq. (j32p is valid. 
Consider next condition (ii), 

UJbf'jl^b) _ _ 1 ^ ^ 

f'{u>b) 7i'/(a;f,) Tra{u!b) 

which is also satisfied when a{uJb) 3> 1; condition (iii) is also obvious from BllIq « TTa{ujb). Finally we find: 

B « — ^e'^ "-'^ " (61) 

The choice of the exponent v is a balance for allowing a maximal range for Ec^. ([32]), which neglects the terms in 
the 3 conditions on equal footing, and the necessity of satisfying the conditions. We expect then v <^ 1. We note that 
the result (|6T|) is closer to the Monte Carlo form^ B ^ jo^ujc than ((59l) above. 



B. Study of the general case 

We present here an analysis of the general case, using the asymptotic expansion in the parameter na. Define the 
function 

F(x) = nrja ^ a„n2e-"'^/'^ (62) 

n 

so that Eq. becomes 

f'iuj)=F[ dwi//(u;i)]. (63) 



The boundary condition for this 1st order equation is /{loc) — ttluc X^n o^nn^ while the condition f'{uJc) — F{0) follows 
from the equation itself. We now generate a 2nd order equation 
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Multiplying by f'{uj) and integrating yields 

-H[f{co)]+H[fico,)]= r f{co,)-^F-^[r{co,)]du;,^-\nlp(. (65) 

Hence in term of the function H{y) — F{x)dx, determined up to one integration constant, one obtains: 

H[f'{u;)] = -HKfiu;)] 

-Hlf'(co^)] 

For the Caldeira-Legget system one can choose H[f'{uj)] = — 7r/'(aj), i.e. a a independent function, which leads to the 
solution in the previous section. In the general case however the function H(y) depends explicitly on a, in the form 
H{y) = ah{y/a) where h is the reciprocal function of F/a. Hence it is not strictly possible to look for a solution of 
the form /(w) = ujg(Kuj) with a a independent g. Explicit integration of (j66p is then required with proper matching 
(j25p at frequency ujq but this will not be attempted here in full generality. For H(y) a power law however, one can 
redefine a scaling function as shown in the next section. 

Instead we will follow an approximate method which is consistent with the one loop RG. The idea is to determine 
the integration constant K by an expansion near ujc where 



This identifies F ^{y) and the function H is then, to 1st order in /'(w) — f'{uJc), 
Using the boundary condition and (|66p 



K = = . (67) 

We can now rederive the RG equation (H^ by a solution of the form f{uj) — ujg{Ku!) such that g{x) does not 
depend explicitly on a. At a; = KlJc we have 



g{x) = ira^Qn 

n 

K{a) 
K'{a) 



I, , K(a) tr-^ 2 /r.c^\ 

^9 yx) -n^a^n (68) 



so that f'{ujc) — V^oiJ2n o-nn^ becomes eq. (|44|) . 

The reasoning below Eq. (j55|) can now be repeated so that /(w) is generated by repeated In embeddings. For the 
effective mass B we need the boundary condition at ujq, i.e. i?[/'(wo)] = H{v'BuJo) = — ln{KBujQ) and g{KuJo) = Buiq 
which yield an equation for the product Bloq, 

gi 75 )=B^o. (69) 

BujQ 

The relation g{Kujo) = Bujq determines then wo and hence, finally, B. Before reaching ujo, at aujo, we expect the 
modification as discussed in the CL system (previous Section), leading to a change in the exponent fi. 

C. Charge Metal system 

In this system we define a mean free path Fermi wavevector kp and theit^i^ the Fourier expansion is identified 

by 

[1 ^ ^] = ^a„ (l - cosnz)/2. (70) 

^4r2sin2(z/2) + 1 
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Hence an ~ ln(r/n) for 1 < n < r, where r = R/l, while a„ w otherwise. Applying <P /dz^ and d'^/dz* at z = 
we get J2n """^ = and J2n ^■n.n'^ = 2r^ + I8r^. 

We show first that the dependence of the effective mass on the radius is B ^ when i? cx), as for a free particle. 
We rely here on the proof of section III, within the variational method, that the effective mass can be found from 
(/) = in Eq. (fTTj) . The action has then the form 

S^nt{0o} = V(l/r)ln(r/n)5{n0o(T)} ^ - / dxlnxS{x9o{T)} 

n=l -^0 

where is a functional of Oo{t), the latter is rescaled as Oo{t) = r6o{T). The action (including the free term Si) is 
then r independent and therefore the effective mass for [^o(''')]'^ is r independent, which after unsealing yields B ^ . 

We proceed to study the variational solution at large r and large a. While realistic metals have a < 1, this study 
supplements MC studies^S,, done at small a. In the region x ^ tt and in the large r limit, the function F{x) takes the 
form: 

F{x) ^ r^F{r^x) , F{y) 2Trr)a dz ln(l/z)z2e-^'^/'^ (71) 



with F(0) — 27rria and F{y) « 2TT^^^riay '^/^ for large y — r^x. For even larger values of x > tt, i.e. y > r^ir the 
function behaves as in the CL regime F{x) — 2£H£££g-2;/7r 

It is useful to define the rescaled function via /(w) = r^/(w) so that in regime x < ir the variational equation 
becomes, for ui < ujc- 

~f{u) = F{r ^) (72) 

which is now independent of r. It can be solved in principle and assuming that the matching frequency cjo occurs in 
this region x < tt we get B = r^B with r- independent conditions /{loq) = Bujq and /'(wq) = v'^^o for 8.iid B. 
Hence luq and B are r independent in the large r limit (they depend on a) and we recover that B ^ . 

Note that in the regime of large y we can use the asymptotic form and the equation (|66p can then be integrated as: 

H[~nuj)] - -27r (27177^)'/' [fioj)]'''' = - ln[if />)] (73) 

where if is an a dependent integration constant. We note that a scaling function can be defined via 

/» = (27rr;a)2/>) = ^9{Ku:) 

where K — K / {2TT'qa)'^ . g{x) satisfies 7r[g(x) + xg'ix)] = ln^[xg{x)], hence g{x) is a independent, except through its 
argument K. 

If ([73|) is used to identify the matching point ([M)) at oja then 



KBUJ^ — ^n'BuJa{2-!Trja)'^/3 

9^T. ^-^e'''^"»(^-''")'/=^) = {27T7^afBu;a (74) 

so that 

{2TnjafBuja = 0(1) (75) 

As we show momentarily, this analysis fails near uja as conditions (i), (ii) fail. As in the CL case, we define uJb 
{u!a < Wfc <C oJc) SO that below ujb the corrected Eq. ([30]) is applied and then we expect luq ~ uJb- 

Alternatively, we can use a scaling form for f{uj) as in Eq. (|57)) i.e. /(w) = 27rwa(w) tJ'{oj) = 27ra;a(a;)r/[a(a;)] . 
Hence (|73)) becomes 

Kuj = ^_e[-'s(-)''('s(-))]^^^ (76) 

which, at a; = Wc, identifies K = e^'^ /{2naujc)- Matching at wq « uJb and using ([7^ 

1 {2TTriaf 



Buq 27ra(wb) 



(77) 
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Note that replacing Wf, — )• Wa recovers Eq. ([75)) . confirming that a(aja) = 0(1). 
For condition (i) we need 



dun TT(2TTrja)'^ 1 . 

- y I ' ^ — (78) 



which is satisfied if a{uJb) ^ 1; clearly at Ua this condition fails. For condition (ii), by a derivative of Eq. (|72p . we 
obtain 

—777 — = o ^"T^^T — ^ F ' < 1 

which is also satisfied when a{uj) <C 1. Finally, condition (iii) is satisfied since B ^ . 
To obtain the effective mass -B, Eq. (|76|) leads to 



where if from (69) is used and k = 27r^/9 at r ^ 1. As in the CL case, we choose a{uj) ~ a'' with 1/ <C 1 so that 
^ Wc, providing a large integration regime for Eq. (32). Finally, the effective mass is 

{2Trria)'^auJc {2Tr)^a^u!c 

It is easy to see that the condition that the frequency belongs to the scale invariant regime and not in the CL 
regime is a: ^ tt i.e.: 

the r.h.s. which can be determined from the solution, depends only on a and not on r, hence this sets a minimum 
radius as a condition. 



VII. CORRELATION FUNCTION 

A. small a perturbation theory 

Independently of the variational method it is also useful to consider the straight small a perturbation theory of 
the action (fTTj) . We consider first the effect of the (/) integration in Eq. ((TT|) . Perturbation expansion in a leads in 
general to a (/> dependence of the form e'^^^<t''^/P where f is a linear combination of the various time variables ti in the 
expansion. The (p integral is then 

/OO 2 2 2 

We expect that the various integrations converge so that at /3 — ?► oo the limit f//? <C if + 0a; can be taken and then 
the sum is dominated hy K — when \(f)x\ < \. We show this explicitely for the 1st order below. 

For 0x = and /3 = oo of (ITT]) we can therefore consider Z^^o- Here we compute the correlation function of cos 9 
order by order in a. The zero-th order is obtained from the free particle action Si ([3]) and given by: 

/-\-oo ^ -j^ 
^(1 - cos(^t))^j^^) = exp(--^M|r|) (82) 

where we have defined ujm = jjj^- To perform the expansion we take the /3 — t- cx) limit in the time integrals since 
these are found to be convergent, while we keep j3 in the dependence. For (px = Q we will rewrite the interaction in 
Eq. (HH): 

Q 1 V- r , , , cos{n[g(r)-g(TO + 27r0(Ti-T2)//3]}-l 
Sint = -^a 2^an I I dTdT — — (83) 



n>l 



oo ./ — oo 
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The first order correction is obtained from the connected average, using Eq. (|8T 

^ 1 

(n - T2) 



cos(0(r)-0(O)))i = ^^a„ / / ^^(r, r_r2)^^^(e'[«M-«(°)+"''(-)-"«(-)l)o,, 



n>l 

where Sk{t,x) = e-f^(.^+('^+^)/l^) /(2MR )^ ^j^g integrals are indeed convergent so that /3 — 00 can 

be taken and SxiTjX) — >■ 1. We have discarded exponentially decaying terms in r such as produced by rt > 1. It 
is important to note that the starting integral is convergent for ti « T2. To see that one can symmetrize in ti,T2 
the term in parenthesis: expansion for n « T2 then yields an additional (ri — r2)^ term. This is a general property 
for all connected averages: the small time apparent singularity is absent. Indeed, expanding the cosine in the vertex 
([55)1 and contracting the two fields with times external to the vertex yields V9{ti)9{t) correlations, which are always 
bounded in the action Si. 

Although integral (|D3[) is tedious to compute, its large t behaviour is easily extracted. It is clear that, assuming 
a Wick decoupling for the e*^ factors, then the exponential decay of each 2-point correlator fixes the value of 6'i — > 
and of 02 — >■ T, leading to a 1/t^ form. To be more precise, the mass term forces the variables ri = + 0(1/ojm) and 
T2 = T + 0(1/wm), i-e. this is the region which dominates the integral in (IDSP at large r ^ 1/ojm- Integration is 
then easy in that region and amounts to replace exp(— iwj\/|ri|) — ^'^^ exp(— iajM|'''2 — t\) — ■^S{t2 — r). 
Note that since the result is only a function of ajj\/T the large r limit is the same as the large ujm limit at fixed r. 
This yields: 

(cos(0(T)-(?(O)))i«aai^^ (84) 

at large r. The amplitude is confirmed by the detailed calculation of the integral given in Appendix [Dl as well as by 
a numerical check. 

We note that this 1 /r^ decay is in general agreement with the constraints derived in Refi^ for the long range XY 
model, very similar to our CL model I. There it was shown that for strictly ferromagnetic LR interactions the spin 
correlation cannot decay slower than the interaction. For the DM model the l/r^ has a ai ^ 1/r coefficient, hence it 
vanishes in the r — ^ c» limit. 



The second order correction can be written as: 



(cos(0(r) - em)2 = ^ E E / - / 

«>1 ri'>l "'^1 "'^4 



^ ^ (85) 



(n - r2)2(r3 - 



At large r the main contribution comes from T2 « r, ra « and ri « T4 (and one deduced by exchange 1, 2 with 3, 4) 
and yields: 

(cos(0(t) - 6{0)))2 « 8a^al^ [ (86) 



'^M Jr, (n 

This integral looks divergent at small times but it is meant to be regularized for ti near zero by the region T2 « r, 
T3 « Ti « T4 « in the above integral ([55]) . This mainly replaces the l/rf factor by a ^ (cos(6'(ti) — 9{0)))i 
factor regular at ri = 0. Similarly, the singularity at ri = r is smoothed by proper integration of (|85p in the region 
T2 w Ti T4 w r, T3 0. Since it is regularized at small times on times of order ~ 1/wa/ the above integral behaves 
as: 

(cos(0(T)-0(O)))2«8a2a?^^ (87) 

To compute the coefficient A we need to perform carefully the integrals in the small time regularization region. The 
question of universality of this amplitude is discussed in Appendix [E) We will not attempt that here but simply note 
that there is no large time divergence in the above integral, i.e. the coefficient A is finite and does not contain any 
log-divergence. 
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B. large q behaviour via matching 

Let us now estimate the correlation function in the large a limit and we restrict to CL for simplicity. For t not too 
large we can just use the straight large a perturbation theory: 

/UJ^ 1 -1-1 
-^(1 - cos(a;T)) — ) « ——^ (88) 
-u:, 27r a\u}\ {uJct)~ 

which presumably is valid only for Iii(uJct) <Si Tra. For larger time one needs to consider renormalization of the 
dissipation. We use the analysis of the previous sections. For large a we expect that we can use the fixed point action 
which is of the form ([3]) with renormalized parameters, i.e. near ac = 1/(0. 742k) and with the mass M replaced by 
the renormalized mass B/R^ with B — B{a). To get an estimate of the correlation function at large r, we can now 
use the above result for the small coupling expansion replacing a by a* and ojm by 1/B[a) (according to Eq. 
(|6ip ). Since a* is not strictly small this will only provide an estimate. One gets, keeping the dominant exponential 
term, 

(cos(^(r) - em) - (89) 
which we expect to be valid for ln(a;cT) ^ n'^a + 0(lnQ:). Eq. (15^)) matches ([55)1 at In sa n'^a. 

VIII. DISCUSSION 

We have studied two types of environments: (i) The Caldeira Legget (CL) system, with relevance to small rings, 
r < 1, or to the Coulomb box problem, and (ii) the dirty metal environment, that can couple to either a charge (CM 
system) or an electric dipole (DM system), with relevance to experiments on cold Rydberg atoms- 

For the CL system, the variational method was shown to be equivalent to an RG scheme, reproducing the known 
two loop result^. Our method provides an expansion to all orders in 1 /a and leads to the renormalized mass Eq. (|61l) , 
which is close to the result of the boundary field theoryi^ and the MC datai^. 

At small a we find a regular expansion without In divergences up to second order, i.e. the RG /3-function for a 
seems to vanish perturbatively. Since we know that the RG flow of a at large a is towards smaller values of a, there 
seems to be three main possible scenarios: (i) the flow towards small a becomes much slower, either exponentially 
due to some putative non perturbative corrections, or to some higher order in a (ii) there is a line of fixed points for 
a < ac with some termination point ac (iii) there is a infinite set of fixed points at small alpha with accumulation at 
zero^. 

In fact the result {cos[9{t)] cos[0(O)]) ^ 1/t^ is a robust one, relating to a theorem on an XY model on a lattice^^. 
That this result is derived in first order in a is remarkable. For large a one should use the scaling to small a and 
then use the former result. 

For the dirty metal problem we show that at large r the whole action scales with r^, leading an r independent 
effective mass B/ E? . Furthermore, we find a scaling form for large r and large a that leads to the renormalized mass, 
Eq. (|79p. For the CM system a < 1 from Eq. ([S]), yet for the DM system a large a may be realized in Eq. (|10l) if 
the dipole has p > e£, i.e. the extension of the Rydberg atom needs to be > £. The large a solution is useful also as 
a complement to the small a MC data at a = 0.19, showing saturation of the effective mass with r. Therefore, the 
claims for an r dependent mass^i are in contrast with both weak and strong a results. We note also that the result 
to first order (cos[6'(r)] cos[6'(0)]) ~ ai/r^ vanishes as ai ^ ^ — suggesting that the nonlinearities associated with 
a become weaker in the large r limit. 

We believe that the correspondence of the variational method with the scaling forms is a useful and instructive 
guide for studies of large variety of nonlinear systems. 

We thank I. S. Burmistrov, A. Golub, P. Guinea, V. Kagalovsky, A. D. Zaikin and G. Zarand for stimulating 
discussions. BH acknowledges kind hospitality and financial support from LPTENS and PLD from Ben Gurion 
University. This research was supported by THE ISRAEL SCIENCE FOUNDATION (grant No. 1078/07) and by 
the ANR grant 09-BLAN-0097-01/2. 
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Appendix A: The rj parameter 

We solve here Eq. (|44| for r]{a). We change variable to x = 1/a and then to y{x) — rj{x)/x 

.'(^) = -TrfiT + — (Ai) 

K[rj[x) — \\ X 

''(^) = <^^' 

therefore 

xy{x) -11, 

- K n — = -n-\ " ^ (A^) 

a;'(y) + = - (A4) 

y 

A general solution of the homogenous part is x = Cie""^ while for a solution to the full equation substitute x = 
A{y)e~'^y so that A'{y) — ne'^y /y, hence 

A{y) ^ K [ —dy = KEi(K2/) + C2 (A5) 
x{y) = (Ci+C2)e-"^ + KEi(K?/)e-'"2/ (A6) 

where Ei is the exponential integral function, with the asymptotic expansion 

Ei(z) = — [l + i + 4 + -] z^oo (A7) 
The boundary condition gives at x 

(A8) 
(A9) 
(AlO) 

where C — Ci + C2- Comparison with Eq. ^ shows that C = 0, a remarkable result. The solution is then 

1 = KaEi{Kar])e-^°''' (All) 
This result is plotted in Fig. 1 for 77 as function of l/{Ka). 

Appendix B: The Log integral 

We present here the mathematical result, i.e. solving the Log integral by using RG and deriving an asymptotic 
solution. Consider the equation: 

f'{y) = -HKf{y)] (Bl) 

TT 

with the boundary condition /(I) — na. This is the equation for the CL system () in the text, defining f{uj) = 
ijj cj iy I '■iJ ci ■ The constant K is parameterized as 

A = e '/Tra 

so that /'(I) — nar]. In general we expect that any pair (0,77) will produce a solution f{y). 
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We rewrite the solution in the form 



f{y) = yg y (B2) 

so that the boundary condition at ?/ = 1 is 

g[k{a)] =Tra. (B3) 

Imagine now a varying boundary condition a and that the function g does not depend exphcitly on a; this imphes 
that r]{a) must be chosen in a specific way. Eq. (|B1[) can be written as 

g{x) + xg'{x) = - Hxg{x)] (B4) 

TT 

where a; = ify. Taking a derivative of (|B3[) yields {dK /da)g' (K) = tt so that Eq. (|B4[) at a; = /^T yields 

a + ---^-^ = ^ ln(A'(a)7ra) = ar](a) (B5) 

leading to a differential equation for r]{a) 

r]-l+ ^^^^ - 1 + (B6) 
a-R''(Q!) TT^ary — 1 + TT^Qf^^ 

Integrating this equation from any initial values (a, 77) yields a function 77(0;); the full solution at j/ < 1 is then obtained 
f{y) = T:ya{y) where the function a{y) is determined by the choice K[a{y)] = K{a)y. Indeed one then has: 

f{y) = ygiKy) = yg[K{a{y))] = TTya{y). 

Therefore one needs to invert the algebraic relation K[a{y)] — K{a)y to find /(y), leading to a form like ([SS)) . In 
general, however, (jB6| is is not easier than the original (jBip . except for the initial values (a = 00,77 = 1) which allow 
an asymptotic expansion with large a. That our physical system satisfies this special tuning is most remarkable. 



Appendix C: integration of RG 

It is instructive to study the one loop RG of the dirty metal system, and compare with the variational solution. 
Consider then the RG equations^, which can also be read off from Eq. ((2T|) 



with have defined an — aa„. Let us consider the function F{x) of Eq. ((62)) setting 7; = 1 there. It becomes now £ 
dependent with: 



dF,{x) ^ F^jx) 
d£ F,(0) 

Change to the variable 



(C2) 



d£ ^^l 

so that 



J^^O)^ = I- (C3) 



dx 

which has the solution 



[^-^]F,{x)^Q (C4) 



F,^^)(x) = F,(o){x + (C5) 
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and with fj, = di/T:Fe{0) we have the general solution 

Pii^) - Poi f (C6) 
Jo 

Note that it has some formal similarity to the variational Eq. if we define £ = —hnj/uc- We proceed to study 
the dirty metal case, with F( = ^^(0), 



Hence by differenciating 



Therefore 



2 2^ SF^ ^ ^ 

6 r p'^i'^ 



2{^-^f'^[Flil-Fl'^]=l = H^) (C9) 



where is a renormalized cutoff. RG terminates at = 1 ^ Fq = 2nar^ , hence 

2(^^^)2/3(2W)i/3^1n(^). (CIO) 

Since powers of r cancel lo^ ~ w^, i.e. the frequency at which the RG is stopped is independent of r, a conclusion 
also obtained in the text. 

Appendix D: calculation of an integral 

The integral given in the text, upon rescaling ajj\/T r is computed as: 

Iff e~5'^l'^l + l'^l~'^^l) -M/ ^ 

' ' .(g^(k-n| + |T2|-|Ti|-|T-T2|) _ J^-J 



-r/2 /-O 



dx^[(l-a;-r/2)(e"-l)-T]+ / da; ^[e-2^(l + r/2 + x) + x - t/2 - 1 

-oo 2; J-r/2 2; 

(■-r/2 -X F-oa -x 

+ dx—[il- e-2-)(l - t/2 + x) - 2a;] + / dx —[(1 - e-^)(l - a; + r/2) - r]} 

Jo ^ Jt/2 X 



/2 

(" _ e3-/2(3^ „ 4)Ei f-^J) - e^/^ (e^V + r + 4) Ei f-^) (Dl) 



2 



-e 

1 



^T2(-21og(T)-27 + l + log(4)) + 0(i3) , r<l (D2) 
^4 + ^ + 0(,e-/2) , ,»i (D3) 

Appendix E: structure of higher orders in small a perturbation 

The discussion of the first and second order corrections in the text suggested that the large time behaviour of the 
integrals could be obtained from a Wick theorem on the e*^ fields with a 6 function correlator, e.g. the structure 
of the second order correction (j85p at large time is an integral dominated by t2 ~ r, ra « and ri « T4, i.e. by 
the region such that the charges in e*"^^'^') should compensate. This was found to generically lead to l/r^ decay. A 
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question is then whether the amphtude of this decay can be obtained to all orders using this property, and how does 
it depend on the details of the short time cutoff. 

Let us first consider a toy model with the same property. It is a gaussian theory of partition sum / Dz{T)e~^ (for 

^ = ^ / ^{r)g{T - T')z*{r') (El) 



where z(r) plays the role of e'^*^"^^ and g{T) « 1/t^ at large t. Being a gaussian variable it reproduces, for a = the 
propagator < z{t)z*{t') >o= 4cl'^/i5(t — r'). There is thus some similarities in the perturbation expansion in a with 
the original model. Here however it is immediate to obtain: 

<z,z;,>=G{t-t') , G(w) = (E2) 

bJM — ^ag{ijo) 

Let us consider two examples for the short time cutoff function in (|E1[) . (i) g{Lij) = 7rajj\/e~'"'/"*^ corresponds to 
a Lorentzian (/(r) = l/(aj^/ + r^). The coefficient of G{t) ~ A/t"^ at large r is obtained from the expansion 
G{uj) = G(0) - ■kA\uj\ + .. and reads A = 8a/((l - 2a7r)a;^) by expanding (|E2]) (ii) ^(w) = -7r|w|e~'^'/'^M which gives 
A = ia/uj\j, i.e. only a first order contribution, all higher orders being zero. 

The above example shows that the amplitude A can depend on the short time cutoff beyond leading order. In that 
case it was however easily calculable by convolutions. To check whether one can indeed predict the 1/t^ coefficient 
more generally, let us consider now the following general discrete XY model of partition sum: 

setting gkk = for convenience. A calculation using mathematica then gives, for x 7^ (here x,p,q belong to an 
arbitrary lattice): 

< e'(''=="^«^ go^ + ^ gop9px + ^ gop9pq9qx - ^gxogL (E4) 

p p,q^O,x 

+ E gxpgpqgqrgrx — ^ gxogoxgopgpx — {gopgxpglx + 9pogpxglp + glxdxpgop) + 0{g^) (E5) 

pz^q^T^p^^^x P7^0,a; P7^0,a: 

Taking as an example a one dimensional chain with discrete t values, and — g{T — r') ^ 1/(t — t')^ at large 
T — t', one sees that up to order 0(g'^) (included) the structure captured by model (jEip is correct to predict the 
coefficient of the l/r^ decay of < e'*-^^"^"^ >. Indeed, up to that order, the l/r^ decay can be obtained from the 
convolution of the g kernel. However some new terms arise at order g* which are not of the above form and do 
contribute to the l/r^ decay, and the calculation of A becomes more complicated than in model (jEl[) . 
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